An Atom Laser is not monochromatic 
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We study both numerically and analytically the possibility of using an adiabatic passage control 
method to construct a Mach-Zehnder interferometer (MZI) for Bose-Einstein condensates (BECs) in 
the time domain, in exact one-to-one correspondence with the traditional optical MZI that involves 
two beam splitters and two mirrors. The interference fringes one obtains from such a minimum- 
disturbance set up clearly demonstrates that, fundamentally, an atom laser is not monochromatic 
due to interatomic interactions. We also consider how the amount of entanglement in the system 
correlates to the interference fringes. 
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Atomic Bose-Einstein condensates (BECs) are often 
referred to as the "atom laser," the matter-wave equiv- 
alent of a laser. Unhke massless photons from a laser 
which do not interact, atoms from an atom laser do in- 
teract with each other, and so far it has not been clearly 
established how and to what extent the monochromatic- 
ity of an atom laser is compromised by such interactions. 
The term "atom laser" is not yet defined universally other 
than in a broad general sense, although there have been 
attempts to be more specific about what an atom laser 
should be, including the directionality of the beam in a 
close analogy with a laser 1]. However, given fundamen- 
tal differences such as the fact that the atoms in a BEC 
are intrinsically stationary, here we shall use the term 
to denote an intense source of coherent atoms for atom 
optical transformations such as refiection, refraction and 
beam splitting. 

In general, the non-zero mass of the atoms necessitates 
rather elaborate and experimentally demanding schemes 
to coherently accelerate and guide the matter-wave. In 
order to mitigate this problem, we propose the use of 
atom optics in the time-domain based on adiabatic pas- 
sage. An important reason for this is that time domain 
atom optics is less prone to decoherence due to, say, atom 
losses as compared to techniques such as guiding a BEC 
in space. Further, the adiabatic passage method is a well- 
known technique in quantum mechanics that provides a 
way to "delicately" control a quantum system via slow 
passage along the energy landscape, while preserving the 
quantum system in one of the energy eigenstates, typi- 
cally the ground state. By minimizing any disturbance, 
and maintaining the condensate in the ground state, one 
can draw conclusions about the intrinsic properties of an 
atom laser unmodified by such additional actions as bal- 
listic expansion. 

To address the question of the nionochromaticity of an 
atom laser, we model both numerically and analytically 
an ideal Mach-Zehnder Interferometer (MZI) for BECs 
in the time domain and deduce the spectral composition 
of an atom laser from the interference pattern generated. 
The construction of a MZI that uses BEC as the source of 
coherent atoms is one of the goals in atom optics owing to 
the promise of, for instance, high precision matter-wave 
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FIG. 1: (a) Schematic diagram of the transformations that 
the potential undergoes over time, starting with beam split- 
ting, addition of a potential step, recombination and final 
beam splitting, (b) Spatio-temporal probability density ob- 
tained by solving the Gross-Pitaevkii equation describing the 
condensate passing through the time-domain MZI. 

based metrology. In a typical optical MZI scheme, co- 
herent light from a laser is passed through a beam split- 
ter which sends the light beam into two paths; one of 
the beams passes through a phase shifter and the two 
beams are then recombined by reflecting off mirrors onto 
a second beam splitter. The phase shift experienced in 
one arm is reflected in the intensity variation at the two 
detectors. We have simulated the time-domain MZI by 
adding an optical potential to the usual quadratic trap- 
ping potential and adiabatically varying its amplitude to 
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FIG. 2: Interference fringes obtained by solving the Gross- 
Pitaevskii Equation (GPE) to model the time-domain MZI 
for various values of the nonlinearity constant g (a) g = 
(b) g = 0.5 (c) (? = 5 (d) g = 10 The sinusoidal fringe for 
the linear g = case is clearly modified for g > i.e. in the 
presence of interatomic interactions, showing that the atom 
laser is no longer monochromatic. 

simulate exactly the action of the optical elements in an 
usual MZI, namely the two beam splitters, mirrors and 
wave guides. The phase shift is provided by adding an 
additional potential step to one "arm" of the interferom- 
eter. Given the adiabatic nature of the method, the two 
arms do not interfere unless they are explicitly recom- 
bined by way of the second beam splitter. The population 
difference AP between the two arms is finally measured 
as a function of the step height. A schematic diagram 
of the required operations is provided in Fig. [TJa). To 
the best of our knowledge, such a scheme has not yet 
been experimentally implemented, although there have 
been several experiments on interference of BEG in op- 
tical potentials Q and in atom chips [1, 0, @ , with the 
important difference that all of these release the conden- 
sates to expand ballistically after the first beam splitter. 
The time-domain MZI proposed here keeps the conden- 
sate in the ground state, better revealing the nature of 
the atom laser. 

The Gross-Pitaevskii Equation (GPE) was used to sim- 
ulate the dynamics of the condensate in the presence of 
the spatio-temporal potential. The GPE has been ex- 
tremely successful in describing the dynamics of BEGs 
under a variety of experimental conditions including op- 
tical latticesjy]. Specifically, we consider the 1-D GPE: 
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where ip{x,t) denotes the condensate mean field 
and V{x, t) denotes the adiabatically changing time- 
dependent potential. As shown later in our analysis, 
the dimensionality of the GPE simulation does not affect 
the final results. We always retain the parabolic trap- 
ping potential, so V{x,t) = ^muJix"^ + Af{t)cos{kx) + 



Vstep{x,Tif4 < t < 2T1/4) where k = 2tt/\, A is the am- 
plitude of the optical lattice potential and T1/4 denotes 
1/4 of the total duration of the simulation T. f{t) pro- 
vides the time dependence for the adiabatically changing 
amplitude of the optical lattice: 
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< t < Ti/4, first beam splitter; 
Ti/4 < t < 2T1/4, first beam splitter; 
2T1/4 < t < 3Ti/4, recombination; 
3Ti/4 < t < 4Ti/4, final beam splitter. 

(2) 

The potential step was introduced during T1/4 < t < 
2T1/4 using a super-Gaussian in one arm Vstep{x,t) = 
/ie-[(*-^i/4)/'^tl'"e-[(^-A/2)/o-,]3"^ ^j^jj ^j^g spatial and 
temporal width of ax — A/3 and at = Ti/4/2 respectively. 
The probability density of the condensate obtained from 
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A = 15 in harmonic oscillator units of the initial trap- 
ping potential is shown in Fig. [Ifb) as a spatio-temporal 
plot. As is clear from the figure, the simulation pro- 
vides exact one-to-one correspondence to the MZI, and 
the adiabatic condition ensures that the wave function 
evolves extremely cleanly. One of the advantages of this 
scheme of interferometry is that the area enclosed by the 
two arms of the interferometer is readily controllable by 
changing A; in particular, with a larger enclosed area, 
higher sensitivity is possible for use in, for instance, in- 
terferometric rotational sensors used in space navigation. 

We plot in Fig. [2] the final intensity difference as a 
function of step height h for different values of nonlin- 
earity g. It is noted that a relatively large phase shift 
is recorded even for small variations in step height, in- 
dicating high sensitivity of the condensate to changes in 
potential height. In the linear case (non BEG) 5 = 
one obtains perfect sinusoidal variation as a function of 
h, while for g > 0, the pattern clearly deviates from the 
sinusoidal, implying that an atom laser loses monochro- 
maticity due to the presence of interatomic interactions. 
This is confirmed by a Fourier transform which reveals 
four closely spaced frequency components. 

In order to better understand this result and to ensure 
that the adiabatic assumptions are not violated, a de- 
tailed analysis of this system can be given as follows. 
First, we assume that the adiabatic passage we used 
above preserves the single mode approximation. The 
condensate after beam splitting and the potential step 
(which imparts relative phase of Q) may then be written 
as: 

il){x, t, 9) = 4>l{x.. d)aL{t) + (t>R{x, d) exp(z6)afl(t) (3) 

where ai(/j)(<:) are bosonic annihilation operators for 
the modes of the matter wave field in the left and 
right arms and 4>l(r){x, d) are the corresponding spatial 
mode functions. These may be approximated by spa- 
tial functions such as shifted Gaussians 4>i^(^ ji){x, d) = 

(7rcr^)~^/^e~(^*'''^/''°'^ that, under adiabatic evolution, 
do not change shape. The phase shift has been mod- 
eled by the explicit insertion of the exp(iQ) term, which 
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is clearly an approximation to the full GPE simulation 
where the phase shift results from the passage of the con- 
densate over a potential step in one arm. The analysis 
is therefore best suited to the cases with lower values of 
g for which interplay between the potential step and re- 
pulsive atomic interactions do not complicate the phase 
relationship given in Eq. 

Since we use the adiabatic passage method, there are 
no additional net physical effects on the BEG undergo- 
ing the remaining two operations that constitute the MZI 
(recombination and then second beam splitting), other 
than those of the interatomic scattering and the natural 
tunnelling flow between the left and right arms. This is 
due to the Josephson effect and the geometry of the sys- 
tem though the tunnelling would be minimized when the 
arms are well separated. The effective interference fringes 
are then given by the difference in the number of atoms 
in the two arms after the system has evolved a certain 
time T through the MZI. With the field decomposition of 
Eq. ([3]) the bosonic Hamiltonian is given by 
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where AE{t) = hu!exp[—(f(t)/a'^] gives the overlap of 
the two spatial modes and is the time-dependent tun- 
neling energy. This overlap clearly depends on the di- 
mensionality of the modes. However, change in di- 
mensionality merely alters this single scaling parameter 
and not the behavior we describe. This Hamiltonian is 
also well-known as the Hamiltonian that describes two- 
component BECs when "L" and "R" denote two dif- 
ferent species instead of the left and right arms. In- 
troducing the Schwinger angular momentum operators, 
J+(_) = "■L{R)°-R(L) and = ^(aiaL - aj^aj^), the 
Hamiltonian is equivalent to a "Lipkin-Meshkov-Glick" 
type Hamiltonian Q, 
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The time evolution operator corresponding to this non- 
linear Hamiltonian has been studied previously [sf. The 
measured intensity difference is given by (Jz(t)) at time 
T, {Jz{t,Q)) = {E)e'"'^RJzR)e~^"'-'k), where R = 



exp 



-iB 



J4 



is the well-known rotation 



operator and H' is the transformed Hamiltonian diag- 
onal in J,:. H' = RHR^ - iR§-^R^ = AEJ,^ - gj^. 
Writing the quantum state in a general form j^*) ~ 
Cm|'Tj)L|Af — m)fi where N is the total number of 
atoms in the system, one obtains. 
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FIG. 3: The interference fringe given by the fractional num- 
ber difference in the two arms of the interferometer, {Jz)/N, 
obtained using the analytical methods under the adiabatic 
assumption. The figure, which shows remarkable agreement 
with the simulation results obtained above using the GPE, 
again clearly indicates the presence of several frequency com- 
ponents in the atom laser from the way it deviates from the 
sinusoidal pattern. 

where Auj = AE/h and 
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i|i?^|r7i') provided 



with the matrix elements TZm^m' = 
by an analytical expression [9|. 

A natural choice of quantum state for the two com- 
ponent BECs with a relative phase is a coherent spin 
state or the Bloch state with equal number of atoms in 
both arms such that \^) = = Tr/2,4> = In this 
case, the coefficients = (C|i„)i/2e*0-'")e/2J where 
denotes the combination, C," = n\/[{n — m)lml]. 
The resulting interference fringe with (Auj)t chosen to 
be {2n + 2)TT and gr = {2n+ ■^)'^, where n is an integer, 
is plotted in Fig. [3] which shows remarkable agreement 
with the result from the GPE simulation. With n « 127, 
we are simulating the g — 0.5 case in the GPE simula- 
tion above, and the effective magnitude of g of the order 
2^ places the system in the Josephson regime. The am- 
phtudes are different from the GPE simulation simply 
due to our choice of (Jz) to be half the atom number 
difference. 

The results confirm that the coherent spin state is a 
natural state for BECs with relative phase arising from 
the adiabatic passage method. With the current system, 
the initial state prepared at the first beam splitter can 
be controlled so that the number of atoms in each arm 
are not equally distributed. For instance, this can be ac- 
complished by adiabatically deforming the potential wells 
asymmetrically to collect more atoms in one arm and 
then continue adiabatically to model subsequent actions 
such as recombination. In such cases, a coherent spin 
state with 9 ^ 'k/2 can be generated, and consequently 
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FIG. 4: Entanglement parameterized by the von Neumann 
entropy of the input state plotted against the visibility of the 
final inference fringes. Inset: Visibility vs. fractional atom 
number difference between the two arms AP of the input 
state. Solid, dashed, and dotted lines: equivalent nonlinearity 
in GPE of g — 0, 0.25 and 0.5 respectively in both the main 
figure and the inset. 

the amount of entanglement between the left and right 
arms is directly controllable. The state with 9 = tt/2 has 
the largest entanglement while the state with 9 = has 
the least amount of entanglement P. 

Finally, how the presence of interatomic interaction 
modifies the connection between the fringe visibility and 
entanglement is of interest from the point of possible ex- 
perimental determination of parameters such as the in- 
teratomic scattering length using this method. The de- 
gree of entanglement between the left and right arms can 
be quantified using the von Neumann entropy E{t) = 
~ iog,(Af+i) T,L=-j |cmP log2 |c™p , which is clcarly con- 
nected to the interference pattern at time r given by 
{Jz{T,'d)). We plot in Fig. [3] the von Neumann en- 
tropy of the input state against visibility V of the fringes 

V = [Imax - Ira%n)/{Imax + Irrnn) whcre / = \{Jz{t))\'^ 

denotes the intensity of the fringe pattern. The visibility 
of the interference fringes is obviously reduced with an 
asymmetric beam splitting that corresponds to ^ tt/2. 



The different line shapes correspond to different values 
of nonlinearity. For a given value of visibility one obtains 
lower entanglement with higher nonlinearity. This ini- 
tially counterintuitive result can be understood from the 
inset of Fig. [H which provides the relationship between 
the fractional number difference between the two arms of 
the initial input state AP and the visibility of the fringes. 
As expected V ^ as AP ^ 1. For AP < 0.75, V al- 
ready reaches 1 as one starts finding atoms in the other 
arm such that ( Jz) starts to take on negative values as a 
function of Q. With higher nonlinearity one finds V ^ 1 
even for larger values of AP. This may be understood 
as the result of higher entanglement due to collisions. 
With higher entanglement, the "wave-like" nature of the 
wave function is enhanced (changing one arm affects the 
other arm instantaneously) which encourages tunneling 
between the two arms, and subsequently higher visibility. 

In conclusion, we have studied a potentially realizable 
scheme for a MZI using BECs. The resulting intensity 
difference between the two arms as a function of phase 
shift is found to demonstrate that, fundamentally, an 
atom laser is not monochromatic but rather comprises of 
a number of frequency components. The adiabatic pas- 
sage method ensured that the only "alteration" made to 
the condensate is that of coherent splitting and the estab- 
lishment of a relative phase between the two arms which 
naturally puts the trapped condensate onto a coherent 
spin state. The dynamics of the split condensate was then 
given simply by a standard many-body Hamiltonian over 
a finite time interval. The fact that monochromaticity is 
lost even in this idealized model implies that one should 
exercise caution when discussing an atom laser as simply 
a matter-wave equivalent of a laser, since additional ef- 
fects such as decoherence due to atom losses is expected 
to degrade the quality even further. Finally we observed 
that with higher nonlinearity, higher fringe visibility can 
be obtained with very uneven beam splitting owing to 
the enhanced entanglement between the two arms. This 
work suggests that there are many new aspects to be ex- 
plored in the time domain quantum control of BEC using 
adiabatic passage. 
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